FORMATION OF SHOCKS IN 
HIGHER-ORDER NONLINEAR DISPERSION PDES: 
NONUNIQUENESS AND NONEXISTENCE OF ENTROPY 

VICTOR A. GALAKTIONOV 

Abstract. Formation of shocks for the third-order nonlinear dispersion PDE 
(0.1) u t = (uu x ) xx in Ix R + , 

is studied by construction various self-similar solutions exhibiting gradient blow-up in 
finite time, as t — > T~ G (0, oo), with locally bounded final time profiles u(x,T~). These 
are shown to admit infinitely many discontinuous shock-type similarity extensions for 
t > T, all of them satisfying generalized Rankine-Hugoniot's condition at shocks. As 
a result, the principal nonuniqueness of solutions of the Cauchy problem after blow-up 
time is inherited, and, under certain hypothesis, any sufficiently general "entropy-like" 
approach to such problems to detect a unique continuation beyond singularity becomes 
illusive. In other words, any attempt to create "entropy theory" for (|0.1|) along the lines 
of the classic one for scalar conservation laws (developed by Oleinik and Kruzhkov in 
the 1950-60s), such as Euler's equation u t + uu x — in M x M + , or by using other more 
recent successful ideas for hyperbolic systems, is hopeless, regardless the fact that these 
PDEs have a number of similar features of formation of shock and rarefaction waves. 



1. Introduction: nonlinear dispersion PDEs and main directions of study 

We study "micro-structure" of shock-type finite time singularities that can occur in 
higher-order nonlinear dispersion equations (NDEs). To describe key features of formation 
of such single point singularities, it suffices to consider the NDE of the minimal third 
order. Its study eventually leads to quite a pessimistic conclusion concerning uniqueness 
of "entropy" solutions and on nonexistence of any hypothetical entropy. 

1.1. NDEs: nonlinear dispersion equations. In the present paper, we continue our 
study began in [5] [S] of basic aspects of singularity formation and approaches to existence- 
uniqueness-entropy for odd-order nonlinear dispersion (or dispersive) equations. The sim- 
plest canonical model is the third-order quadratic NDE (the NDE-3) 

(1.1) u t = {uu x ) xx =uu xxx + 3u x u xx in Rx(0,T), T > 0. 
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We pose for (ll.ip the Cauchy problem (the CP) with locally bounded and integrable 
initial data 

(1.2) u(x, 0) = u (x) in R. 

It is principal that we consider the CP, where the solution u(x, t) is supposed to be defined 
by initial data (11.21) only, plus, of course a generalized Rankine-Hugoniot-type condition 
on the speed of propagation of shocks, which follows from the equation integrated in a 
shock neighbourhood. In other words, we are assuming that the CP does not require 
any a priori posed conditions on the shock wave lines (though of course such ones exist 
and can be determined a posteriori). Otherwise, with such conditions, we arrive at a 
free boundary problem (an FBP) for the NDE (11.11) . which requires other mathematical 
methods of study and can be well-posed (unlike the CP). We will touch possible FBP 
settings for (II. ip . 

The physical motivation of the NDEs such as (11.11) and other odd-order nonlinear PDEs, 
which appear in many areas of application, with a large number of key references, are 
available in surveys in [7J § 1] or [SJ § 1] , so we omit any discussion on these applied issues. 

1.2. NDEs in general PDE theory. For future convenience, we briefly present a more 
mathematical discussion around NDEs; cf. [8, § 1]. In the framework of general theory 
of nonlinear evolution PDEs of the first order in time, the NDE (II. ip appears the third 
among other canonical evolution quasilinear degenerate equations: 



(1.3) u t = —\ {u 2 ) x (the conservation law), 

(1.4) u t = \ {u 2 ) xx (the porous medium equation), 

(1.5) u t = \ {u 2 ) xxx (the nonlinear dispersion equation), 

(1.6) u t = — | (\u\u) xxxx (the 4 th -order nonlinear diffusion equation). 



In (ll.6p . the quadratic nonlinearity u 2 is replaced by the monotone one \u\u in order to 
keep the parabolicity on solutions of changing sign. The same can be done in the PME 
(II .4p . though this classic version is parabolic on nonnegative solutions, a property that is 
preserved by the Maximum Principle. Further extensions of the list by including 

(1.7) ut = -\ (u 2 ) xxxxx (theNDE-5) and 

(1.8) u t = \ (\u\u) xxxxxx (the 6 th -order nonlinear diffusion equation), etc., 

are not that essential. These PDEs belong to the same families as (11.51) and (11.61) respec- 
tively with similar covering mathematical concepts (but more difficult in some details). 

Mathematical theory of the first two equations, (ll.3p (see detailed survey and refer- 
ences below) and (ll.4p (for quoting main achievements of PME theory developed in the 
1950-80s, see e.g., P Ch. 2]), was essentially completed in the twentieth century. It is 
curious that looking more difficult the fourth-order nonlinear diffusion equation (11.61) has 
a monotone operator in H~ 2 , so the Cauchy problem admits a unique weak solution as 
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follows from classic theory of monotone operators; see Lions [T3l Ch. 2]. Of course, some 
other qualitative properties of solutions of (11 .6p are more difficult and remain open still. 

It turns out that, rather surprisingly, the third order NDE ( 11. 5ft is the only one in this 
basic list that has rather obscure understanding and lacking of a reliable mathematical 
basis concerning generic singularities, shocks, rarefaction waves, and entropy-like theory. 



1.3. Mathematical preliminaries: analogies with conservation laws, Riemann's 
problems, and earlier results. As a key feature of our analysis, equation (ll.ip inherits 
clear similarities of the behaviour for the first-order conservation laws such as ID Euler's 
equation (same as (11.31) ) from gas dynamics 

(1.9) u t + uu x = in lxR + , 

whose entropy theory was created by Oleinik [TIT [To] and Kruzhkov [12] (equations in 
M. N ) in the 1950-60s; see details on the history, main results, and modern developments 
in the well-known monog raphs UJ El HH0 

As for (11.91) . in view of the full divergence of the equation (11.11) . it is natural to define 
weak solutions. For convenience, we present here a standard definition, mentioning how- 
ever that, in fact, the concept of weak solutions for NDEs even in fully divergent form 
is not entirely consistent and/or very helpful, to say nothing about other non-divergent 
equations admitting no standard weak formulation at all; see [8]. 

Definition [Til. A function u = u(x,t) is a weak solution of (li.ip . (QUI) if 

(i) u 2 G Ll oc (R x (0,T)), 

(ii) u satisfies in the weak sense: for any test function <f(x,t) G C^°(M x (0, T)), 
(1-10) jj uipt = \ jj u 2 tp xxx , 

and (iii) satisfies the initial condition in the sense of distributions, 

(1.11) ess lim^o j u(x, t)ip(x) = j u (x)ip(x) for any ip G C£°(R). 

The assumption T < oo is often essential, since, unlike (11.91) . the NDE (11.11) can produce 
complete blow-up from bounded data, [8j § 4]. Thus, again similar to (jl.9p . one observes 
a typical difficulty: according to Definition [TJl, both discontinuous step-like functions 

/-, io\ o i \ ■ ( 1 for x > 0, 

(1-12) S T (*) = T sign* = T j_ lfor x<Q) 

are weak stationary solutions of ( 11. II) satisfying 
(1.13) (u 2 ) xxx = 



1 First study of discontinuous shocks for quasilinear equations was performed by Riemann in 1858 |18j 
(by Riemann's method); see [21 HZ] for details. The implicit solution of the problem (|1.9|) . u — uo(x-ut) 
(containing the key wave "overturning" effect), was obtained earlier by Poisson in 1808 [IB] : see [T7] . 
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in the weak sense, since in ( 11.101) u 2 (x) = S?L(x) = 1 is C 3 smooth (x — does not count). 
Again referring to entropy theory for conservation laws (11.91) . it is well-known that 

U-(x, t) = S-(x) is an entropy shock wave, and 
u+(x,t) = S+(x) is not an entropy solution. 

This means that 

(1.15) U-(x, t) = S-(x) = — signx 

is the unique entropy solution of the PDE (II. 9p with the same initial data S-(x). On 
the contrary, taking S + initial data yields the rarefaction wave with a simple similarity 
piece-wise linear structure 

{—1 for x < —t, 
~ for |x| < t, 
1 for x > t. 

The questions on justifying the same classification of main two Riemann's problems 
with data (jl,12p for the NDE (11.11) and to construct the corresponding rarefaction wave 
for S + (x), as an analogy of (11.161) for the conservation law were addressed in [H [6]. This 
was done by studying the following self-similar solutions of (II. lft : 

(1.17) u.(x,t) = f(y), y=(zf73 => Uf')" = \f'y in K, ffroo) = ±1. 

Here, by translation, the blow-up time in (I1.17P reduces to T = 0. It was shown that, in 
the sense of distributions or in L\ oc , 

(1.18) u-(x,t) -»• S-(x) as *->0". 

Therefore, S^(x) is a 5-entropy shock wave (see concepts in [6J), while S + (x) is not and 
creates a typical rarefaction wave given by the global similarity solution 

(1.19) u + (x,t)=F(y), y=-^ (FF')" = — ~ F'y in R, F( T oo) = T l, 
where indeed F(y) = —f(y), so that, in the same sense as in ()1.18p . 

u + (x, t) — > S+(x) as t — > + . 

1.4. Layout of the paper: new shock patterns and nonuniqueness. In Section [21 
we study new shock patterns, which are induced by other similarity solutions: 

(1.20) u4x,t) = (-t) a f(y), y=j3p, P = where a G (0, |) and 

(1.21) (//')" -(3fy + af = in R_, /(0) = /" (0) = 0. 

The anti-symmetry conditions in (11.211) allow to extend the solution to the positive semi- 
axis {y > 0} by —f(—y) to get a global pattern. The case a < in (11.201) . corresponding 
to the strong complete blow-up was studied in detail in (SJ § 4] in the parameter range 

(1.22) a G [ — ^j,0), where a c = — ^ is a special critical exponent. 

For convenience, we revise some of these blow-up results in the range (ll.22p in Section [2j 
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Obviously, the solutions (11.171) . which are suitable for Riemann's problems, correspond 
to the simple case a = in (ll.20p . We prove that, using positive a, allow to get first 
gradient blow-up at x = as t — > 0~, as a weak discontinuity, where the final time profile 
remains locally bounded and continuous: 



1.23) u_{x,0' 



Cq\x\ p for x < 0, 

— 

-Cq\x\ 13 for x > 0, 

where Co > is an arbitrary constant. Note that § < 1 for a < |. 

Therefore, the wave braking ("overturning") begins at t — 0, and, in Section [31 we show 
that it is performed again in a self-similar manner and is described by similarity solutions 

(1.24) u+(x,t)=1fF(y), y = $, (3 = ±±2 where 



(1.25) 



(FF')" + (3F'y-aF = in R_, 

F(0) = F >0, =CoMf(l + o(l)) as y -> -oo, 



where the constant Co > is fixed by blow-up data (11.231) . The asymptotic behaviour 
as ?/ — > — oo in (11.251) guarantees the continuity of the global discontinuous pattern (with 
F(—y) = —F(y)) at the singularity blow-up instant t = 0, so that 

(1.26) u-{x,0~) = u + {x,0 + ) in R. 

Rather surprisingly, for a fixed Co > in (I1.23P obtained by blow-up evolution as 
t — > _ , we find infinitely many solutions of the extending problem (11.25p . This family of 
solutions $cb = {F(y; F ), F > 0} is a one-dimensional curve parameterized by arbitrary 

(1.27) F = F(0) > 0, 

which does not have any clear boundary or extremal points. We also show that for F = 
(11.251) does not have a solution. In other words, this solution family does not contain any 
"minimal", "maximal", or "extremal" points in any reasonable sense, which might play a 
role of a unique "entropy" one chosen by introducing a hypothetical entropy inequalities, 
conditions, or otherwise. 

A first immediate consequence of our similarity blow-up/extension analysis is as follows: 

(1.28) 



in the CP, formation of shocks for the NDE (II. ip leads to nonuniqueness. 



The second conclusion is more subtle and is based on the mentioned above fact on the 
homogeneous structure of the functional set <3>cv if ^c — { u +( x it), F £ ^c } is the 
whole set of weak solutions of (11.11) with initial data (ll.23p . then, for the Cauchy problem 

for onp, 



1.29) there exists no general "entropy mechanism" to choose a unique solution. 



Of course, we cannot exclude a hypothetical situation, when there exists another, non- 
similarity solution of this problem, which thus does not belong to the family *f>c = { u +} 
and is the right solution with a proper entropy-type specification. In our opinion, this is 
suspicious, so that we claim that (I1.28P and (11.291) show that the problem of uniqueness 
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of weak solutions for the NDEs such as (11.11) cannot be solved in principal. On the other 
hand, in a FBP setting by adding an extra suitable condition on shock lines, the problem 
might be well-posed with a unique solution, though proofs can be very difficult. 

In other words, non-uniqueness in the CP is a non-removable issue of PDE theory for 
higher-order degenerate nonlinear odd-order equations (and possibly not only those). In 
fact, the non- uniqueness of solutions of (11.251) has a pure and elementary dimensional 
nature. Indeed, this is a third-order ODE with the general solution depending on three 
parameters, which are more than and excessively enough to shoot the right behaviour 
at infinity and hence to get the continuity (11.261) at the blow-up time. This guarantees 
non-uniqueness of the extension for t > and is a core of this difficulty. 

Increasing the order of the PDE under consideration, we then enlarge the dimension of 
the parametric space, and this surely will imply even stronger non-uniqueness conclusions. 
For example, non-uniqueness and non-entropy features are available for the fifth-order 
nonlinear dispersion equation (the NDE-5) 

(1.30) u t = —{uu x ) xxxx in Ix M + 

and others; see basic models in [6j [7J . 



2. Gradient blow-up similarity solutions 

In this section, we consider the blow-up ODE problem (11.211) . Actually, this ODE is 
not that difficult for application of standard shooting methods, which in greater detail 
are explained in [HI § 3,4]. Moreover, even similar fifth-order ODEs associated with the 
shocks for the NDE-5 (11.301) also admit similar shooting analysis [7J, though, in view of 
the essential growth of the dimension of the phase space, some more delicate issues on, 
say, uniqueness of certain orbits, become very difficult or even remain open. Therefore, 
in what follows, we will omit or even will not mention some technical details concerning 
the problem (ll.2ip and recommend consulting |SJ in case of troubles. We will widely use 
numerical methods for illustrating and even justifying some of our conclusions. For the 
third-order equations such as (11.11) . this and further numerical constructions are performed 
by the MatLab with the standard ode45 solver therein; see more details in [H § 3,4]. 

Let us begin with a simpler fact concerning this problem and some simple asymptotics 
for matching purposes. We recall the elementary symmetry of the ODE (11.211) 

(2.1) //^-/, 

which allows us to put two conditions at the origin. Such solutions have a sufficiently 
regular asymptotic expansion near the origin: for any A < 0, there exists a unique solution 
of the ODE (OB . satisfying 

(2.2) / ( y)=A y + i^ 3 + i^iy 5 + ... as y - 0. 

The uniqueness of such asymptotics is traced out by using Banach's Contraction Principle 
applied to the equivalent integral equation in the metric of C(—8, 5), with 5 > small. 
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We use the following scaling invariance of the ODE in (11.211) : if fi(y) is a solution, then 

(2.3) faiy) = a3 fi{a) * s a solution for any a ^ 0. 

By scaling (12.31) . the parameter A < can be reduced to a single value, say A = — 1. 

Let us describe the necessary bundles of the 3D (in fact, 4D, for the non-autonomous 
case) dynamical system (11.211) . Firstly, due to the scaling symmetry (12.31) . there exists 
the explicit solution 

(2.4) f*(y) = my 3<Q for y < ~ 1 - 

The overall bundle about (12.41) is 2D, which is obtained by the linearization: 

(2.5) f(y) = f.(y)+Y(y) [My)y}"'-f3Y'y + aY = ... = 0. 
The linearly independent solutions are 

(2.6) Y(y) = y m (m - 3)(m 2 + 9m - 20a - 2) = 0. 

Hence there exist two negative roots m± composing a 2D stable manifold. Fortunately, 
f*(y) is strictly negative for y < and attracts no orbits from the positive quarter-plane. 
Secondly, the necessary behaviour at infinity of f(y) is: 

(2.7) f(y) = C \y\^ (l + o(l)) as y - -oo (ft = J), 

where Cq > is a constant, which can be arbitrarily changed by scaling (12.31) . It is 
important for the future conclusions to derive the whole 3D bundle of solutions satisfying 
(12. 7p . As usual, this is done by the linearization as y — > — oo: 

f(y) = fo(y)+Y(y), where f (y) = C (-y)7 

=► C [(-y)fYY" + {3Y'(-y) + aY + l[f (y)y" + ... = 0. 

According to classic WKBJ-type asymptotic techniques in ODE theory, we look for solu- 
tions of (12.81) in the exponential form with the following characteristic equation: 



(2.9) Y(y) ~ e<~^\ 1 = l + \{l-f ) >l =^ C ( 7 a) 2 = -ft a± - . 



; 1 
7 V CO" 



This gives the whole 3D bundle of the orbits (12. 7ft : as y — > — oo, 
(2.10) /(y) = C (-y)f + (-y) 5 [C lS in (I J^ G {-yV) + C 2 cos (I J^yT)] + 



where Ci j2 G R. The slow decaying factor (— y) 5 in the double scale asymptotics (12.101) is 
not essential in what follows, so we do not specify the exponent 5 < therein. 

The behaviour ( 12.91) and (12.101) give crucial for us asymptotics: due to ( 12. 71) . we have 
the gradient blow-up behaviour at a single point: for any fixed x < 0, as t —>■ - , where 
y = x/{-t)P -> -oo, 



(2.11) = (-t) Q /(y) = (-t) a Co|^|^(l + o(l)) ^ C |x 

uniformly on compact subsets, as required by (11.231) . 
The geometry of the above asymptotic bundles yields: 



3a 
1+a 
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Proposition 2.1. For any fixed A < in (12. 2 j) . the problem ( \1.21h admits the unique 
shock wave profile f(y), which is an odd function and is strictly positive for y < 0. 

Uniqueness follows from the asymptotics (12.21) and scaling invariance ( 12. 3ft . Global 
existence as infinite extension of the unique solution from y = 0~ follows from the structure 
and dimension of the bundles (12. lOj) (it is 3D, i.e., comprises the whole phase space of the 
equation in this quadrant). The 2D "bad" bundle, composed from the orbits (12.51) . is not 
connected with the anti-symmetric ones. Concerning positivity, which is rather technical, 
see some details in [HI § 3,4]. 

Figure [T] shows a general view of similarity profiles f(y) for various values of the pa- 
rameter a G [0.01, 0.5] with the fixed derivative 

/' (0) = -10. 

In Figure El we show enlarged oscillatory structure corresponding to (12.101) of some of the 
profiles closer to the origin. Figure [3] continues to explain the oscillatory behaviour (12.101) 
of solutions on different y-scales. 

The next Figure H] shows how the oscillatory features dramatically increase for negative 
a. For a = —0.099, which is very close to the critical value —j^ in (ll.22p . we observe a 
"saw-type" profile of maximally allowed oscillatory structure; see more details in |8, § 4]. 

Finally, we claim that, besides odd blow-up similarity profiles, there exist other not that 
symmetric, for which the conditions at the origin in (ll.2ip do not apply. This construction 
can be performed similar to that in [H § 3,4]. The blow-up solutions constructed above 
are sufficient for our main purposes. 
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Blow-up similarity profiles for a =0.5, 0.4, 0.3, 0.2, 0.1 , 0.05, 0.01 




FIGURE 2. Enlarged oscillatory behaviour of f(y) from Figure [T] for y G [—80,0]. 



Blow-up similarity profiles for a = 0.05, 0.01 , and 



Blow-up similarity profiles for a = 0.05, 0.01 , and 





( a ) V 



-500 



(b) y 



-1000 



FIGURE 3. Enlarged oscillatory behaviour of f(y) from Figure Q] for y € 
[-500,0] (a) and y G [-1000,0] (b). 



2.1. On self-similar collapse of shocks. It is curious that the similarity solutions 
( 11.201) given by the ODE (11.211) can describe collapse as t — > 0~ of shocks. These are 
given by profiles f(y), which instead of the anti-symmetry conditions in (jl.2ip . satisfy 

(2.12) /(0) = / >0, / / (0) = / 1 <0, and /"(0) = f 2 e R. 



Since the bundle at infinity (12.101) is 3D and hence exhausts all the trajectories there, 
shooting with the parameters (12. 12)) yields an orbit f(y; f , fx, fa), which, for a wide 
range of /o,i,2, has the behaviour (12. 7p with a Cq > and creates the data (12. lip as 
* -> 0". 
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Blow-up similarity profiles for a = 0.01 , 0,-0.05, and -0.099 
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y 

FIGURE 4. Oscillatory behaviour of f(y) for positive and negative a. 

In addition, for the corresponding similarity solution (ll.20p . the shock disappears since 

(2.13) [u-{0,t)) = 2f (-t) a -f as t -+ 0~ . 

However, since the data at t = 0~ has a typical form (12.111) . a new shock will be created 
as t — > + , which we are going to explain in Section [3j where the parameterization such 
as in (12.121) for the ODE (11.251) will be used in greater detail. 

2.2. On smooth rarefaction waves. These are global in time solutions of the form 
(11.241) . with the ODE as in (ll.25p . Then taking F(y) = —f(y), one observes how the 
weakly singular initial data 

(2.14) -u-(x,0~) given in §TM) 

collapse into the smooth (even analytic) solution u + (x,t) for t > 0. This is quite similar 
to the same phenomena as in [S] and [B], and we will not comment on this anymore, but 
indeed will use the global solutions (ll.24p . however in their discontinuous hypostasis. 

3. NONUNIQUE SIMILARITY EXTENSIONS BEYOND BLOW-UP 

3.1. Nonuniqueness of similarity solutions. As we mentioned, a discontinuous shock 
wave extension of blow-up solutions (11.201) . (11.211) are performed by using the global 
ones (ll.24p . (11.251) . Actually, this leads to watching a whole three-parametric family of 
solutions parameterized by their Cauchy values at the origin: 

(3.1) F(0) = F > 0, F'(0) = Fi < 0, and F"(0) = F 2 eR. 

The 3D phase space for the ODE in (11.251) has two clear stable "bad" bundles: 
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(I) positive solutions with "singular extinction" in finite y, where F(y) — > as y — > 
y$ < 0. Indeed, this is an unavoidable singularity following from the degeneracy of the 
equations with the principal term FF'" leading to the singular potential ~ ^; 

(II) positive solutions with the fast growth about the explicit solution: 
(3.2) F*(y) = -> +oo as y -> -oo. 

The 2D stable bundle is similar to that in (12.51) . 

Both sets of such solutions are open by the standard continuous dependence of solutions 
of ODEs on parameters. The desired solutions are situated in between those two stable 
open bundles; cf. arguments in P, § 3,4]. 

As usual, it is key to derive the whole 2D bundle of solutions satisfying (12.71) . Similar 
to (12.81) . we perform the standard linearization as y — > — oo in (11.251) : 

f(y) = F (y) + Y(y), where F (y) = C (-y)% 

C [(-y)7Y]"'-(3Y'(- y )-aY + l[F (y)}'" + ... = 0. 

The WKBJ method now leads to a different characteristic equation: 

(3.4) Y{y) ~ e<~yy, 7 = 1 + \ (l - f) > 1 =}► C ( 7 a) 2 = /3,a ± = ±± 
so that the only admissible root is a_ < 0. This gives a 2D bundle of the orbits (12. 7ft : 

(3.5) f{y) = C {-y)% + (- y yC 2 e a -(-y^ + ... as y - -oo, del. 
As a result, we have the following: 

Proposition 3.1. Let a G (0, |). For any fixed F > and Fi < in (13.11) . t/iere exzsfe 
a unique F 2 G R sucn t/iat t/ie problem l \ 1.25b has a solution F*(y) for some Co > 0. 

The proof is performed by shooting as in [SI E] by using the stable bundles indicated in 
(i) and (ii) above. Such techniques are currently well-established for various higher-order 
ODEs. As a similar and more complicated example of a fourth-order ODE, we refer to 
the methods in [IT] , where by shooting technique existence and uniqueness of a positive 
solution of the radial bi-harmonic equations with source: 

(3.6) A 2 r u = u p for r = \x\ > 0, u(0) = 1, u'(0) = u"'(0) = 0, u(oo) = 0, 
was proved in the supercritical Sobolev range 

P > Psob = fzf, where N > A. 

Here, analogously, there exists a single shooting parameter being the second derivative 
at the origin F 2 = u"(0); the value Fq = u(0) = 1 is fixed by obvious scaling. Proving 
uniqueness of such a solution in [11] is not easy and lead to essential technicalities, which 
the attentive Reader can consult in case of necessity. Note that, instead of the global 
behaviour such as (I3.2p . the equation (13. 6p admits the blow-up one governed by the 
principal operator 

u {A) + ... = u p (u -> +oo). 
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Shooting similarity profile F.(y) fora=0.3: F'(0)=-3.398..., F(0)=1, F"(0)=10 




-20 -18 -16 -14 -12 -10 -8 -6 -4 -2 



FIGURE 5. Shooting a proper solution F*(y) of (ll,25h for a = 0.3 and data 
F(0) = F = 1, F"(0) = F 2 = 10, with Fi = F'(0) = A* = -3.398.... 

The solutions vanishing at finite point otherwise can be treated as in the family (I). 

Thus, as a result, we obtain a two-parametric family of solutions of fl 1 . 2 5 [) with an 
arbitrary fixed C > parameterized by F > and F x < 0. For a fixed constant C > 
(uniquely given by the blowing up limit t — > the family is one-parametric: 

(3.7) $c = {F«(y;F Q ), F > 0}, 

which for convenience we parameterize by the value Fq = F*(0) that measures the jump 
of the shock. The actual parameterization of the family (13 .7p is not of importance. For 
instance, in Figure O we show shooting F*(y) for a = 0.3, which actually explains the 
strategy of proving of Proposition 13.11 In Figure [6l the same is done for a = 0.1 and 
smaller -F(O) = 0.1, F"(0) = 0. In Figured again for a = 0.1, we show a couple of more 
shooting with respect to parameters F 2 (a) and Fi in (b). 

3.2. Remark: on nonexistence of solutions with Fq = 0. This is also a principal 
issue. Indeed, if a solution of (11.25P for F = 0, i.e., F(y) having the asymptotics near 
the origin similar to (12. 2p . 

(3.8) F{ y) = A y-^y 3 + ^^\y 5 + ... as y^0 (A < 0), 

would exist, then the similarity solution (11.241) would describe a smooth collapse of the 
initial singularity (ll.23p . and would actually mean an extra nonuniqueness in the problem. 
Fortunately, this is not the case and the single parameter A < in (I3.8P (actually reducing 
to A = —1 by scaling (12. 3p ) is not sufficient to shoot the necessary asymptotics as y — > — oo 
given in the second line in (ll.24p . i.e., the corresponding asymptotic bundles as y — ► 0~ and 
y — > — oo are not overlapping. We do not prove this carefully (but indeed this can be done, 
since the "non-overlapping of bundles" is large enough), and, as a key illustration, present 
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Shooting profile F.(y) for a=0.1 : F'(0)=-0.555480988..., F(0)=0.1 , F"(0)=0 




Ql 1 1 1 1 

-1 -8 -6 -4 -2 



y 

FIGURE 6. Shooting a proper solution of (|1.25j) for a = 0.1 and data 

F(0) = F = 0.1, F"(0) =F 2 = 0, with F 1 = F'(0) = A* = -0.55548098.... 

Shooting similarity profile F.(y) for a=0. 1 : F'(0)=1. 13285..., F(0)=1, F'(0)=-1 Shooting similarity profile F.(y) fora=0.1: F'(0)=-0.257714..., F(0)=1 , F"(0)=0 




(a) F = l,Fi = -1 (b) F = 1, F 2 = 



FIGURE 7. Shooting a proper solution F*(y) of (jl.25p for a = 0.1 and data 
F(0) = F Q = 1, F'(0) = Fx = -1, F"(0) = F 2 = A* = 1.13285... (a) and 
F(0) =F Q = 1, F"(0) = F 2 = 0, F'(0) = Fi = A* = -0.257714... (b). 

Figured! where this nonexistence is carefully (with Tolerances about 1CU 10 ) checked for 
a — 0.1 (a) and a = 0.3 (b). 

However, for more complicated nonlinear PDEs, such as, e.g., [7] 

v>u (uUxjxxxxi etc., 

where the resulting ODEs are higher order with multi-dimensional phase space, by some 
reasons, such solutions can be available, meaning that formation of shocks can be re- 
versible and is not a unique option. 
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Nonexistence of "regular" F (y) for a=0.3: F(0) e [-6,-1 ] 



Nonexistence of "regular" F (y) for a=0.1 : F(0) e [-6,-10 ] 





(a) a = 0.3 



(b) a = 0.1 



FIGURE 8. Towards nonexistence of a solution of (|1.25p . (|3.8j) . with Fq = 0, for 
a = 0.3 (a) and a = 0.1 (b). 

3.3. Further discussion around nonuniqueness and non-entropy issues. Thus, we 
have explained the required non-uniqueness H 1 . 2 8 1) of the solution of the Cauchy problem 
(11. ip . (11.231) : taking any profile G $c yields the self-similar continuation (ll.24p . with 
the following behaviour of the jump at x = 0: 



(3.9) 



-[u + (x,t)}\ x=0 = -(u + (0 + ,t)-u + (0-,t)) = 2F t a >0 for t > 0. 



Rephrasing the result, let us emphasize that, in the similarity (i.e., an ODE) representa- 
tion, it has a pure dimensional origin: the problem (11.251) has too many solutions. More 
precisely, for the given 3D shooting via parameters {F Q , F\, F 2 } in the 3D phase space of 
the ODE in (ll.25p . the dimension ID of the non-suitable bundles in (I) and (II) is not 
sufficient to define a unique solution F*(y), up to scaling (12.31) as usual. In other words, 
the desired uniqueness could be achieved if 

(3.10) the "bad" bundles in (I) and (II) are 2D in the 3D phase space 

(plus certain natural structural "transversality" as a genetic property). Hence, (I3.10p is 
a geometric recognition of a possible uniqueness (and entropy that assumes extra hard 
work) extension of singular shock wave solutions. 

Note that since these shocks are stationary, the corresponding Rankine-Hugoniot (the 
R-H) condition on the speed A of the shock propagation: 

(3.11) A __K^ 1 ' 



\x=0 







is valid by anti-symmetry. As usual, this condition is obtained by integration of the 
equation (11.11) in a small neighbourhood of the shock. Alternatively, (13.111) is obtained by 
approximating the solution via a travelling wave (T W) 



(3.12) u(x,t) = f(x-Xt) 



(//')" 
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-A/ = (//')' =► A 



[(//')'] I 
If] ly=o' 



which coincides with ( 13. lip . Recall again that the R-H condition does not assume any 
novelty and is a corollary of integrating the PDE about the line of discontinuity. 

Moreover, the R-H condition ( 13. lip also indicates another origin of nonuniqueness: a 
symmetry breaking. The point is that the solution for t > is not obliged to be an odd 
function of x, so we can define the self similar solution (11.241) for x < and x > using 
different triples of parameters {F^, F^, F^}, and the only extra condition one needs is 
the R-H one: 

(3.13) A = =► [(FF')']=0, i.e., + (F^f = F+F+ + (F+) 2 , 
which can admit other solutions than the obvious anti-symmetric one 

F =-F+, Ff = F+, and F 2 =-F+. 

Overall, in any case, we have at least a one-parameter family of solutions of the CP 
tyc = { u +i x it)} f° r t > depending on arbitrary parameter Fq > 0. This family 
is not discrete and hence does not reveal any particular solution u + , exhibiting some 
special properties as being, say, maximal, minimal, or extremal in any sense. Under a 
natural assumption, this confirms the extra negative statement fl 1 . 2 9 [) on nonexistence of 
a sufficiently general entropy-like inequality, condition, and/or a procedure to detect a 
unique solution, as a "better" special one. At least, if even an "entropy-like" procedure 
would have been derived (somehow, by any hypothetical means) , in view of the nonunique 
formation of shocks everywhere, the resulting solution would have never been relative to 
data in any metric, i.e., Hadamard's well-posedness concept would have been violated 
anyway. 

Indeed, different regularization of the NDE can lead to different solutions, i.e., by using 
a parabolic regularization: 

(3.14) u e : u t = (uu x ) xx - eu xxxx (e > 0), 

with the same data Uq(x). Proving that the regularized sequence {u e , e > 0} is a compact 
subset in some metrics is indeed a very difficult problem, which remains open for general 
data (see [5] for some details concerning such problems). Anyway and however, assuming 
that a suitable compactness of {u £ } is already available, we do not think that a proper 
unique solution u(x,t) can be obtained as the unconditional limit 

(3.15) u(x, t) = lim e _^ + u e (x, t). 

More precisely, in view of the above nonuniqueness of post-blow-up similarity extensions, 
we believe that (13.151) has infinitely many partial limits and, along various subsequences 
{^k} + , the corresponding sequences {u £k } can converge to different solutions u + with 
various values of rescaled shocks F > 0. We justify this in an example below. 

Example. We take a = | and Co = 1, so that the blow-up initial data (11.231) are 

(3.16) u (x) = y\x\ for x < 0, 

with the odd extension for x > 0. Performing in (13.141) a natural scaling 

(3.17) u £ (x,t) = eh{y,r), y = J^, r = -4^ 
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deletes e and yields the following uniformly parabolic problem: 



(3.18) V T = (VVy)yy ~ V yyyy frl IR_XlR + , V Q ( Xj ) = \J \ JJ \ frl R_, 

with the anti-symmetry conditions v = v yy — at y — 0. It is not a great deal, by using 
classic parabolic theory, to prove global existence of a unique classical solution of the 
CP (13.181) . though it takes some efforts. But the main issue is not existence/uniqueness. 
Indeed, according to (I3.17p . the study of the convergence (I3.15P as e — > + reduces to 
delicate asymptotic behaviour of v(y, r) simultaneously as r — > +oo and y — > — oo, 
which represents a very difficult problem. However, the above detected nonuniqueness 
of the similarity extensions makes such an open asymptotic problem excessive and not 
that necessary (quite fortunately it seems). Note that studying the behaviour of the 
hypothetical proper limit (13.151) at the point of discontinuity x = assumes detecting 
the limit of {u £ (xk,t)} (t > small is fixed) for rather arbitrary sequences {xk} — > 0~ as 
£ — > 0. We then claim that such a limit essentially depends of the choice of {sk} — > + , 
i.e., the limit of the sequence 

i 

(3.19) u £k (x k ,t) =elv(^,^) as k -> oo 

£ k £ k 

is very much {xk}- and {efc}-dependent. 

We also expect that using other types of regularization, e.g., by classic Bubnov-Galerkin 
methods of finite-dimensional approximations on suitable functional Riesz bases (see e.g., 
strong applications in Lions' classic monograph p3]), which obviously give globally ex- 
isting solutions, will also lead to similar nonuniqueness issues in the limitH, so these are 
unavoidable difficulties of modern PDE theory. 

In other words, we then expect that, for a number of higher-ordeiH nonlinear PDEs 
with singularities such as blow-up, extinction, or shock wave formation, fully consistent 
(see above) general entropy-like procedures to reveal a unique solution cannot be available 
in principle. As we have shown, this is just prohibited by a sufficient dimension of the 
phase space (depending on the spatial order > 3 of the PDE), which allows to shoot a 
continuous family of solutions beyond singularity, whose family does not have any isolated 
and/or boundary points. 

3.4. On uniqueness for the FBP setting. Evidently, the uniqueness can be restored if 
an extra condition at the shocks is assumed, which poses an FBP for (11. ip . For instance, 
following Figure [6] and 0(b), this happens if we fix 

(3.20) F 2 = F"(0)=0 =^ u 5KB (0 ± J t) = 0. 

2 Therefore, the issue of the Galerkin uniqueness, which can be treated as a simplest idea for an 
entropy-like construction, becomes author's unrealizable dream. 

3 This is a principle issue: for general second-order nonlinear parabolic equations with blow-up and 
obeying the Maximum Principle, there always exists a unique proper minimal solution, which does not 
depend on the type of monotone regularization of the equation; see [TU1 § 2] or [U Ch. 6]. 
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Indeed, the uniqueness (in the present self-similar setting) is restored if the set (13.71) 
contains a unique such profile, i.e., 

(3.21) $ Co n {F;'(0) = 0} = {F m (y), with a fixed F > 0}. 

Then this fixes the unique "shock divergence" (13.91) with F = F beyond the singularity. 
In a PDE setting, existence-uniqueness of a solution of the FBP (11.11) . (I3.20p is a difficult 
problem, which, for some simple geometric configuration of shocks, can be solved by 
traditional FBP methods, such as von Mises transformations and others. In general, 
mathematical difficulties can be extremely challenging. 

A more general free-boundary condition on shock lines can be predicted from the struc- 
ture (ll.24p of (generic) similarity continuation beyond blow-up. One should only take 
into account that parameters a and (3 cannot enter such conditions, since these essen- 
tially depend on a priori unknown blow-up "initial data" u(x, 0~). Since at the shock at 
x = 0, by 

(3.22) u + =t a F , (u + ) x =t a -?F u and (u + ) xx = t a ~ 2 ^F 2 , 

it is easy to reconstruct a general FBP condition at shocks, which suits for arbitrary a: 

(3.23) uu xx = k(u x ) 2 , where sGl, 

Since both sides are of order t 2a ~ 2/3 , this condition well corresponds to any of similarity 
formation of shocks (what happens for other configurations, is another delicate story). 
Then (I3.20p is obtained for k = 0. Formally, k = oo yields Fi = 0, i.e., a kind of 

a "Neumann" FBP : u x = at shocks. 

Of course, it is necessary to check, for which k the condition 

(3.24) F F 2 = ^F,) 2 

yields a unique profile F* 6 <3>c* , and this occurs for arbitrary a. As customary, posing 
necessary free-boundary conditions is an applied physical issue, though checking and 
predicting the well-posedness of the FBPs occurred is indeed a mathematical problem. 

4. Final remarks: on origin of uniqueness for the ID Euler equation 

In this connection, it is interesting to trace out similar origins of the uniqueness in the 
CP for the Euler equation (11.91) . Its obvious advantage is that it is solved via character- 
istics and the solutions are given by an algebraic relation: 

(4.1) T = if == ' > x - ut = const. =>- u(x , t) = u (x - u(x , t)t) . 
Let us assume the same initial data (11.231) . so that for x < we have the equation 

a 

(4.2) u = Co(ut — x)0 , where now, dimensionally, j3 = 1 + a. 

Obviously, setting here x = 0, corresponding to the permanent position of this stationary 
shock, yields the unique value of the solution at the shock: 

s 

(4.3) u = C {ut)% u{0~,t) = Cg- a t0=^ = C$t a for (3 = 1 + a. 
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Or, analogously, the same uniqueness is guaranteed by the fact that, for the first-order 
ODEs corresponding to rescaling of (ll.9p . the phase space is also one-dimensional, that 
allows a unique (and very simple) matching of two bundles. The corresponding similarity 
solutions for t < and t > respectively are given by: 

u4x,t) = (-trf(y), y = j^ w => ff' + Pf'y-af = 0, 

u + (x,t) = t a F(y), y = % FF' — (3F'y + aF = 0, 

where (3 = 1 + a, and then as in (14.31) we have to have 

(4.5) u_(0~,t) = C(jV for t > 0, i.e., F(0) = C%. 

Of course, the ODEs in (14. 4p are explicitly integrated. As usual, the last one for 
the post-blowing up behaviour, is responsible for a unique continuation beyond blow-up. 
Integrating it yields the unique solution F*(y) defined by the algebraic equation 

(4.6) F> = F = -yP, yP' = =► = C* for y < 0, 

where Cq is the constant in (14. 2 j) (a similar procedure applies to the blow-up profile /). 
It then follows from (14.61) that the rescaled shock jump is also uniquely determined: 

(4.7) F = F*(0) = eg, 

and indeed this precisely gives the time behaviour (14. 3 j) and (14. 5 p of the shock value. 

Note that, unlike calculus in (13.131) . a symmetry breaking is obviously impossible here. 
Indeed, the R-H condition uniquely implies the symmetry: 

(4 R) A - [F2] - (Ill^l = F o+ F o - n =^ F+ - F- 

A -2[F]~ 2(F+-F -) - 2— - U =^ *0 - ^0 • 

It is not an exaggeration to say that this unique micro-structural extension of generic 
blow-up singularities at any point eventually reflects a true dimensional and ODE origin 
of those classic entropy theories constructed for conservation laws in M by Oleinik in 
the 1950s and by Kruzhkov in 1960s for equations in M. N . Indeed, if, for the simplest 
ODE similarity problem for (11.91) . the matching would be nonunique and non-discrete 
(plus something else), these very influential theories would have never been appeared and 
created. We thus claim that this nonuniqueness and non-entropy are precisely the case 
for the Cauchy problem for the NDE (II .ip and can be expected for a number of other 
nonlinear dispersion (and not only those) higher-order equations from PDE theory of the 
twenty first century. 
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Shooting similarity profile F.(y) fora=0.3; F'(0)=A.=-0.94108..., F(0)=1, F"(0)=0 




